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The vector and the axial-vector meson couplings with the vector and the axial-vector currents 
respectively at finite temperature have been obtained in Ref. [1] by calculating all the relevant 
one-loop Feynman graphs with vertices obtained from the effective chiral Lagrangian. On the other 
hand, the same couplings were also derived in Ref. [3] by applying the method of current algebra and 
the hypothesis of partial conservation of axial- vector current (PCAC). The latter method appears 
to miss certain terms; in the case of the vector meson coupling with the vector current, for example, 

\^ , a term containing the pun coupling is missed. A similar situation would also appear for the nucleon 

f"^ ■ coupling with the nucleon current. In this note we resolve these differences. 

o 

The thermal two-point functions of the vector and the axial-vector currents are evaluated in Ref. [1] from their 
. one-loop Feynman graphs with vertices given by the effective chiral Lagrangian of strong interactions [2]. At low, but 
not too small, temperature T, the p meson coupling F p with the vector current, for example, is found to depend on 
temperature as 
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I. INTRODUCTION 



Here F n is the familiar pion decay constant, — 93MeV and gi is the pojir coupling, gi = .87. Although there is no 
strict power counting of momenta according to the number of loops when heavy mesons are present, the above result 
is expected to hold to within an accuracy of about 20% [2] . 

Earlier, the authors of Ref. [3] also evaluated the same quantities in a somewhat different way. At low temperature 
they expanded the thermal trace over states in the two-point functions, retaining only the vacuum and the one pion 
state. The forward pion-current amplitude so obtained was apparently evaluated with the methods of current algebra 
II 1 and PCAC, reducing it to a combination of vacuum two-point functions of the vector and the axial-vector currents. 
Oh! It allowed them to conclude that the coupling Fj would change with temperature according to Eq.(l.l), but without 
the g\ term. 



A similar situation would arise also in the calculation of the temperature dependence of the nucleon coupling A 
with the nucleon current [4]. It has been calculated in Ref. [5] as 



by evaluating in the effective theory the one-loop Feynman graphs for the thermal two-point function of nucleon 
currents. Here g^ is the axial-vector coupling constant of the nucleon, g^ = 1.26. If, however, we follow the procedure 
of Ref. [3] to work out the same quantity by applying apparently the current algebra techniques, we would miss the 
g\ term in Eq.(1.2). 

The evaluation of the thermal part of one-loop Feynman graphs is actually a tree level calculation, as the corre- 
sponding part of the propagator of a particle (here pion) contains a mass-shell delta function. Now, at the tree lavel, 
the current algebra plus PCAC approach is strictly equivalent to the method of chiral perturbation theory [6] . One 
thus wonders why the two approaches would lead to different results for the temperature dependence of the couplings 
stated above. 

Historically, the idea of PCAC was introduced [7] to extract the pion pole contribution to matrix elements of 
operators containing the axial- vector current. At small momentum k carried by the current, this is the only source 
of dominant ciontribution, as far as processes involving only pions (Goldstone bosons) are concerned. However, when 
there are heavy (non-Goldstone) particles in the process, there is another potential source of dominant contribution 
to such matrix elements, namely the graphs, in which the axial-vector current attaches to an external heavy particle 
line. The additional heavy particle propagator so created can be singular at small momentum k. In fact, this is the 
mechanism contributing dominantly to processes of soft pion emission and absorption [8]. 
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In this note we rederive the temperature dependence of Fj and A T in the framework of current algebra in the chiral 
symmetry limit. We begin with the Ward identities for the relevant amplitudes involving the pion and the currents. 
We show that there are indeed heavy particle poles in the above two-point functions that are as dominant as the pion 
pole at small momentum k. In the case of the vector current, it is the co meson pole in the neighbourhood of the p 
meson pole, when we consider temperatures higher than the mass difference of the two mesons. In the case of nucleon 
current, it is the nucleon pole itself, there being no mass difference here to deal with. It is these contributions, giving 
rise to the g\ and the g\ terms in Eqs. (1-1-2) respectively, that are ignored in Ref. [3]. 

In Sees. II and III we present current algebra derivations of Eqs. (1.1) and (1.2) respectively, where the contributions 
of the heavy particle poles are also calculated explicitly along with that of the pion pole. Our discussions are contained 
in Sec.IV. ' 



II. CURRENT ALGEBRA DERIVATION OF EQ.(l.l) 

The derivation is based on the thermal two-point function of vector currents, V*(x),i = 1,2,3, with which the p 
communicates. At low temperature /3 _1 (= T), the heat bath is dominated by pions. Thus the leading correction 
to the p meson propagation is expected to come from its scattering off the pions. Accordingly the thermal trace is 
expanded as [3], 

Tre- fiH TV^(x)Vj(0)/Tre- fiH 

= (0\TV;(x)Vj(0)\0) + [ J2(* m (k)\TV*(x)VJ(0)\* m (k)} conn , (2.1) 

J (Zir) i\k\ m 

where n(k) = (e^' — 1) _1 is the distribution function of pions. T also denotes the time ordering of all operators that 
follow it. (• • -) ccmn is the connected part of the matrix element. The first term gives the p pole in terms of its mass 
and coupling in vacuum. To find their changes at finite temperature to leading order, we use the method of current 
algebra to extract the leading terms of the pion matrix element at small k. 




FIG. 1. Amplitude W^TJ" 3 (q,k,q' ,k') describing the scattering n(k') + V v (q') — > A a (k) + V M (q). Isospin indices are omitted 
in the diagram. 

The appropriate Ward identity follows from the momentum space amplitude 

W^{q,k,q',k') = f d 4 zd 4 xd 4 x'e 4fc - z e^-e-^-\0|TA™(z)V?(x)^(x')k"(fc')) , (2.2) 



describing the scattering process shown in Fig.l. As usual, we contract it with k a and integrate by parts with respect 
to z and ignore the surface term to get 

k a W™™ 3 (q,k,q',k') = i J d 4 zd 4 xd 4 x'e*- z e 4 «- K e - 4 «'- x '^(0|TA™(z)V;(x)V; J (x')k"(fc')) . (2.3) 

Since the axial- vector current is conserved in the chiral limit, the divergence operator gives non-zero contributions 
only when it acts on the theta functions in the definition of the time-ordered product, giving 

d«(0\TA™(z)Vi(x)Vj(x')\n n (k')) 

= S(z x o )(0|T([A"(z), V;(x)]Vj(x'))\n n (k')) + S(z a - x' o ){0\T([A™(z), VJ(x')]V;(x))\n n (k')). 
The equal time commutators of currents are given by current algebra, 

K l (z), V;{x)] ZQ=X0 = te mU Al(x)S 3 (z - x) . 
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Also using translational invariance to remove one more integration in Eq. (2.3), we get the Ward identity as 
= -(2?0 V(9 + k — q' — k') 
[e ml1 J d 4 xe l{q+k) - x (0\TAl{x)Vj(0)\ir n {k')} + e mjl J d 4 xe tq - x (0\TV^x)Al{0)\TT n (k'))} . (2.4) 

To find the leading (constant) piece in k a W™™i in the limit k a — > 0, we need only to calculate the terms in W™™3 
that are singular in this limit. One such term is, of course, the tt pole contribution, as shown in graph (a) of Fig. 2, 
that can be readily evaluated by inserting the pion intermediate state in the matrix element of Eq.(2.2), 

k <* W ™nij -(2*f5\q + k-q'- k')F n J d 4 xe iqx (tt™ (k) \TV* (x)Vj (0) |tt" (*')) , (2-5) 

where F ff is defined by 

(0|^(«)|7r n (A!)) - i5 mn F„k a e- ikz . 

As we shall show below, there is another finite contribution to k a W™™^ as fc" ^ due to the u> pole in the vicinity 
of the p pole. Denoting this contribution as 

fc Q W™ (^) 4 S 4 (q + k-q'- k')k a R™i(q, k, k'), 
the Ward Identity (2.4) gives the result 

J d 4 xe iq ' x (Tr m (k)\TV* (x) Vj (0)|7r n (fc')> 

= e md J d 4 xe l{q+k) - x (0\TA l fl (x)VJ(0)\Tr n (k')) + e mjl J d 4 xe iq - x {0\TV^(x)Al(0)\w n (k')) + k a R™™ j • (2.6) 

The first two matrix elements on the right of Eq.(2.6) can be reduced to vacuum matrix elements, again by using 
Ward identities. Thus for the first one we define the amplitude 



wj^(q,k,k')= / d 4 z'd 4 xe 



ik' -z i{k+q)-x 



(0\TA l ^x)V>$)AW)\0) 



(2.7) 



One proceeds exactly as before by contracting w™™ 3 with k' 13 . Here the pion pole is the only dominant contribution 
at small k! , so that the Ward identity for this amplitude becomes 



F v J d 4 xe^ k+q > x (0\TAl(x)VJ (0)\TT n (k r )) 

= e nl1 ' J d 4 xe iq '- x (0\TV^(x)Vj(0)\0) + e n]l ' J d 4 xe l{q+k) - x (0|T A l ^{x) A l l {0)\Q) . 



(2.8) 



A similar equation holds for the second matrix clement on the right of Eq.(2.6). We thus get from Eq.(2.6) the pion 
matrix element in Eq.(2.1) as 



d 4 xe iqx ]T (iT m (k)\TV* (x) Vj (0) | tt™ (fc)> 

m 

= -Jll d4 ^ 9 ' X < l T ^W(0)|0) + -^y d 4 xe* q - x (0\TAl(x)Al(0)\0) + jrk a R^(q,k) . (2.9) 




n I V v 1 ~ / \ K-' ^ It- 

fa) (bi) ' <b 2 ) 

FIG. 2. Singular contributions to W™^ 3 (q, k, q' , k') as k, k' — ► 0, when the vector currents are coupled to p meson. 
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The remaining task is to evaluate the u> meson pole contribution, shown in graphs (61) and (b 2 ) of Fig. 2. However, 
we do not evaluate them as Feynman graphs with the intermediate lines as propagators, but insert intermediate states 
in the spirit of the method of current algebra, as we already did with the pion pole extraction in Eq.(2.5). Noting the 
quantum numbers of u, it is given by 

T, I fe«|TiC(»)^(*)M*.^»M«.^)W WKW>3T3 + "O"" 1 lerm ' <2 ' 10) 

Its correction to the p pole is obtained by isolating this pole at the vector current vertices. Then extracting the 
coordinate dependence by translation invariance, we can integrate over the coordinates to get delta functions in 
momentum, which are then removed by intergals over the intermediate momenta. We thus get 

^ E (0|^(0)|p b (g,a 2 ))(p"( <Z , ( T 2 )| J 4™(0)|c( g + fc,a 1 )> 

<Tl,CT2,C3,6,C 

(w(q + k, a,)\^{k'), p c (q>, a 3 )) (p c (q\ a 3 )\Vj(0)\) ^ , 2 £ - fc)2 _ ^ } + crossed term . (2.11) 

To evaluate the matrix elements in Eq.(2.11), we write the relevant pieces of the effective Lagrangian based on the 
chiral symmetry of QCD. In the notation of Ref. [2], these are 

1 /V <Wr> + ^iW<^{^^}}- (2-12) 



2V2 m v + 4 

Here the fields are all 3 x 3 octet matrices. denotes the vector meson fields and V^ u their field strengths. The field 
strengths (and f^) are those of external fields, v 1 ^ and a* , with which the vector and the axial- vector currents are 
coupled in the QCD Lagrangian. u M is related to the covariant derivative of the pseudoscalar (Goldstonc) fields. All 
these matrix valued fields transform only under the unbroken subgroup of QCD. Taking the trace denoted by (• • •), the 
above terms are actually invariant under the full chiral symmetry group. Restricting to terms needed in our calculation, 
Eq.(2.12) gives immediately the vector and the axial-vector currents in the effective theory [ud^v = {d p u)v — ud^v], 

V; = F p m pP ^ , Al = -gre^A^dV* , 

and also the pujir interaction, 

i-pum = tt^aXci p O TT , 

where the symbols denote the corresponding fields. Then we have for the matrix elements 
(0|T/;(0)|p b (<z, <7 2 )) = S ib F p m p e p (q, a 2 ) , 

(p b (q, a 2 )\A™(0)\u(q + k, aj) = i6 mb 9l e^ Xa (2q + kYe x *(q, a 2 )e»(q + k, u\) , 

where ^(q, a) is the polarisation state of the vector meson with momentum q and spin projection a = (1, 0, —1). The 
remaining matrix element is a S-matrix element, which to first order in perturbation expansion is, 

(w(« + Mi)k n (*0.P c (9>3)> 

= (2tt) 4 % + k-q'- k')5 nc i^e^ Xa (2q + 2k - k') v k'° e X (q' \ <T 3 )e»* (q + k, a{) . 

r-rt 

Inserting these matrix elements in Eq.(2.11), we can carry out the sums over the polarisation states of the vector 
mesons. After a little simplification, we get 



p k a R^ v 4 ^(j?) (5 2_^2 )2 { (5 + A . ) 2_ m 2 + iq _ k) 2_ m 2}' (2-13) 



where the Lorentz tensor is given by 
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A^((7,fc) = e aPlt _ t ef j , Yu q l3 q l3 'k~ ( k~ ( ' 

= (q-k) 2 g^ + q 2 k li k u - (q.k)(k^q v + k v q^) , 
in the chiral limit (k 2 = 0). Wc can now write the complete result of our calculation as 

i J d 4 xe iq - x Tre- ^ 'TV*{x)Vj {0)/Tre~ 0H 

= { 1 ~%) l j d4 ^(0|T^(x)^(0)|0) + ^ J dW'-*<0|T4(s)Ai(0)|0) 
9i Y (q 2 -ml) f d 3 kn(k) A^^fc) 



V (q 2 -mj) ffk 
'■) (q 2 ~m 2 +ie) 2 J (2tt) 



F*J {q 2 -m 2 p + ie) 2 J (2n) 3 2\k\(q 2 -ml + te) 2 -4(q- k) 2 ' 
where we reinstate ie to define the poles and 

d 3 kn(k) T 2 



(2.14) 



J 



r d 6 h 

~ J (2tt 



(27r) 3 2|fc| 24 



To proceed with the evaluation, we have to take note of the kinematics of the thermal two-point functions. Any 
such function can be decomposed into two invariant amplitudes. We now set q = with q n = E, when the two 
amplitudes as functions of E become related to each other. Thus it suffices for us to consider, for example, the trace 
of the two-point function in its Lorentz indices. The free p pole term in vacuum then reads as 



•; J d 4 xe iqx (0\TV*(x)VM(0)\0) 3 ^ 2 ^ p ^ )2 . (2.15) 



Next we find this trace for the third term on the right in Eq.(2.14) in the vicinity of the p pole. If we ignore the 
mass difference between p and u> mesons, it can be immediately evaluated to reproduce the g\ term in Eq.(l.l). In 
the realistic case with Am = m w — m p = 12 MeV [9], it is useful to note that the fc-integral is essentially cut off at 
|fc| = T, because of the pion distribution function present in it. Thus if we are interested in the behaviour of FT for 
temperatures larger than Am, as is usually the case, we can again ignore Am to leading order to get Eq.(l.l). But 
for temperatures small compared to Am, the fc-integral is 0(T A ) and the g\ term will not appear in Eq.(l.l). 



III. CURRENT ALGEBRA DERIVATION OF EQ.(1.2) 



The properties of the nucleon can be studied by constructing a nucleon current r]%(x) out of the quark fields, having 
the quantum numbers of the nucleon [4]. Here a, b, ■ ■ ■ and A, B, ■ ■ ■ denote respectively the isospin and the Dirac 
spinor indices. To study its properties at low temperature, we again expand, as in the case of vector currents, the 
thermal two-point function of nucleon currents, 

Tre-? H T V a A (x)f] h B (0)/Tre-P H , (3.1) 

in terms of the vacuum and the one pion state, and apply the method of current algebra to estimate the resulting 
pion matrix clement of these currents, 

^(^"(fc)|T^(x)^(0)|7r m (fc)), (3.2) 

m 

for small pion momenta in the chiral limit. As expected the calculation is similar to the earlier one with vector 
currents, except for the fermionic nature of the nucleon current. 
The required Ward identity follows from the amplitude 

Za,ABiP,k,p',k')= J d 4 zd 4 xd 4 x'e lk - z e lp - x e" lp - x \0\TA^(z)r]%(x)f] b B (0)\TT n (k')} (3.3) 

Again using the axial-vector current conservation and the equal time commutation relation, 

[A" \z), V a A (x)] Zo=Xn = ( l5 ) AcV l,(x)6 3 (z-x), 
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we follow the steps as before to get the Ward identity, 
k a Z™7B = -*(2n) 4 5\p + k-p'-k') 

m \ac / m \ cb 

—J (j 5 )ac J d 4 ^+ fe )-(0|T^(x)^(0)|^"(fc')) + J d A xe^ x {Q\T lA (x)f] c c {Q)W\k')) [-) ( l5 ) CB 



(3.4) 



The dominant contribution to Z™ 7 }^ at sm aU k, k' arise from the pion and the nucleon poles. Denote the latter 
contribution as 

ZZTb (2n) 4 S\p + k-p'- k')SZ n A a B(P, k,pf) . 
The forward scattering amplitude (3.2) can now be obtained as 

d 4 xe^ ]T ( fc ) \ Tr )A {x)n B (0) k m (k)) 

m 

= / d 4 xe^(0\T VA (x)f] h B \0) " 2^2/ rf 4 ^-(0|T^(x)^|0) + jrk a S a a b AB (p,k) , (3.5) 

where fj — 7577. 

It remains to evaluate the contribution of the nucleon intermediate state to Z™^. Of course, it must be inserted 
twice more to couple it to the nucleon currents. We get 

ZZTb^ E (0\v a A (0)\N d (p,a 2 ))(N d (p,a 2 )\A^(0)\N c (p + k,a 1 )). 

f 1 jf 2 jf 3 ,c,d,e 

(AT C(j3+ fc;CTl )l^(y^3),^(^)>(^(^^3)l^(0)|0) (p2 _ ^ )(p , 2 _ + fc)2 _ ^ } + c r o S scd tor m (3.6) 

The coupling of 77 to nucleon is given by the matrix element 

(0\r lA (0)\N d (p,a))=S ad Xu A (p,a), 
where u(p, a) is the Dirac spinor of the nucleon. The piece in the effective Lagrangian 

gives the nucleon contribution to the axial-vector current and the nNN coupling. We then have 
(N d (p,a 2 )\A™(0)\N c (p + k,cn)) = \g A u{p,a 2 ) la ^u{p + k,cj 1 )x d] T m x c 



(N c (p + k, a^N^p', a 3 ),ir n (k')) = (2n) i 5 4 (p + k-p'-k') (~^r) 0(p + k, <JiW^u(p', a 3 ) X ct r n X e , 

where x's denote isospin states of the nucleon. Inserting these matrix elaments in Eq.(3.6) and carrying out the spin 
and isospin summations, we get 

STab = -^[0 + ™h°0 + t ™) W + m)] AB 

i 3 

^ x ^ T m T n b + crossed term (3.7) 

(p 2 - m 2 N )(p /2 - m 2 N ){(p + k) 2 - m 2 N } 
We evaluate it in the forward direction for p = and set po = E to get 

±k"S a a b = -6S ab (^A\. k {i> + rnm + rn) *' „ f _— ^— — + __L__") (3.8) 



in the vicinity of the pole. Inserting the results (3.5) and (3.9) into an equation similar to Eq. (2.1), the corrections 
change the free pole term in vacuum 

to one with A replaced by A T given by eq.(1.2). 
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IV. DISCUSSION 



In this work we use the method of current algebra to derive the temperature dependence of certain couplings, that 
have already been worked out with chiral perturbation theory [1,5]. Basic to the former approach is the extraction of 
terms of the appropriate amplitudes containing the axial- vector current operator, that are singular as the momentum 
carried by this operator goes to zero. Generally speaking, the use of PCAC can incorporate in the calculation only 
a part of such singular pieces, namely those due to the pion pole. But there can be other singular pieces in the 
amplitude from heavy particle poles also, that must be calculated separately. The calculation of of the authors of 
Ref. [3] is incomplete in that it does not take the latter contribution into account, explaining the discrepancy between 
their results based apparently on current algebra and PCAC, and that of the effective theory. 

In the case of the two-point function of vector currents, the contribution from the ui meson pole can be as dominant 
as the pion pole in the neighbourhood of the p meson pole. To be specific, it brings in a small scale Am in the 
expression. For temperatures small compared to Am, the u> meson contribution is 0(T 4 ), changing the mass and 
residue only to this order. As T increases past a narrow transition region around T = Am, the behaviour changes 
over to <3(T 2 ), when Fj is given by Eq. (1.1). A similar result holds also for the axial-vector meson a\ [1]. For the 
two point function of nuclcon currents, it is the nucleon pole that must be included. Here the situation is simpler, as 
there is no mass difference in the problem. Accordingly A T is given by Eq.(1.2) for all temperatures. 

It is interesting to compare here this old method of current algebra to its modern development, namely that of 
effective field theory. The Ward identities used in the former method relate ampltudcs with different number of pions, 
while in the effective theory one calculates the amplitudes directly in perturbation expansion. When one is in need 
of vertices involving heavy particles in a current algebra calculation, he has, however, to turn to the effective chiral 
Lagrangian. The insertion of intermediate states is rather clumsy when compared with the use of propagators in the 
effective field theory. Also the two terms in Eqs. (2.13) and (3.8) can be obtained neatly from the mass shell delta 
function in the temperature dependent part of the pion propagator [1,5]. Thus even at the tree level, at which current 
algebra works, one cannot but appreciate the technical and conceptual superiority of the method of effective field 
theory. 
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